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ABSTRACT 

We numerically investigate the effect of feedback from the ionizing radiation heating from massive 
stars on the evolution of giant molecular clouds (GMCs) and their star formation efficiency (SFE), 
which we treat as a time-dependent quantity. We consider the GMCs' evolution from their formation 
by colliding warm neutral medium (WNM) streams up to advanced star-forming stages. We find, 
in agreement with our previous studies, that the star- forming regions ("clouds") within the GMCs 
are invariably formed by gravitational contraction, so that their internal non-thermal motions must 
contain a significant component of global convergence, which does not oppose gravity. After an initial 
period of contraction, the collapsing clouds begin forming stars, whose feedback evaporates part of the 
clouds' mass, opposing the continuing accretion from the infalling gas. The competition of accretion 
against dense gas consumption by star formation (SF) and evaporation by the feedback, regulates 
the clouds' mass and energy balance, as well as their SFE. We find that, in the presence of feedback, 
the clouds attain levels of the SFE that are consistent at all times with observational determinations 
for regions of comparable SF rates (SFRs). However, we observe that the dense gas mass is larger 
in general in the presence of feedback, and that the total (dense gas + stars) is nearly insensitive to 
the presence of feedback, suggesting that the total mass is determined by the accretion, while the 
feedback inhibits mainly the conversion of dense gas to stars, because it acts directly to reheat and 
disperse the gas that is directly on its way to forming stars. 

We find that the factor by which the SFE is reduced upon the inclusion of feedback is a decreasing 
function of the cloud's mass, for clouds of size ~ 10 pc. This naturally explains the larger observed 
SFEs of massive-star forming regions. We also find that the clouds may attain a pseudo-virializcd 
state, with a value of the virial mass very similar to the actual cloud mass. However, this state differs 
from true virialization in that the clouds are the center of a large-scale collapse, continuously accreting 
mass, rather than being equilibrium entities. Finally, we also calculate the density probability density 
functions of the clouds, finding that they in general exhibit the bimodal shape characteristic of ther- 
mally bistable flows, rather than a lognormal form, which is characteristic of isothermal flows. This 
supports suggestions that low density, atomic gas pervades molecular clouds. We conclude that the 
general state of star-forming regions is likely to be one of gravitational collapse, although most of the 
mass in a GMC may be no t participating of the instantaneous star- forming activity, as recently sug- 
gested bv lElmegreenl (|2007n : that is, that SF is a spatially and temporally intermittent phenomenon, 
with strong, localized bursts interspersed within much more quiescent gas. 
Subject headings: interstellar matter - stars: formation - turbulence 



1. INTRODUCTION 

The evolution of giant molecular clouds (GMCs) is a 
key ingredient in our understanding of the star formation 
process. In particular, the low observed star formation 
efficiency (SFE) at the scale of whole GMCs (~ 2%; 
iMvers et al.lll986l ) remains a topic of strong debate, with 
there being two main competing scenarios that attempt 
to explain it. These scenarios refer essentially to the 
effect of the stellar feedback (mainly from massive stars) 
on the star-forming clouds. One is the scenario that the 
stars quickly disrupt their parent clouds by dispersal 
and/or photoionization, before the gaseous mass of the 
cloud is completely converted to stars (e.g. iWhitworthl 
19791: lElmegreenl fl983l: iFranco. Shore fc Tenorio-Taglel 
1991 IWilliams fc McKed 119971 : 

Hartmann. Ballestcros- Paredes fe Berginl 1200 ID . In 



this case, the SFR of active star-forming sites may be 
large for brief periods, and then halted by the very stars 
that have just been formed. 

In the other scenario, the role of stellar feedback 
is to drive turbulent motions within the GMC, which 
oppose its self-gravity, allowing it to remain in near 
hydrostatic equilibrium for times significantly l onger 
than its free-fall time (tfr) (iKrumhqlz fc McKed 120051: 
Krumholz. Matzner. fc McKed 120061 : iLi fc Nakamura 
2006). In this case, the low efficiency of star formation 
would be due to the dual role of supersonic turbulence in 
self-gravitating clouds, of opposing global collapse of the 
cloud while promoting local collapse of turbulent density 
enhancements, whi ch involve small fractions o f the 
total cloud mass dKlessen. Heitsch fc Mac Low! 120001 : 
Vazquez- Scmadcni, Ballestcros -Paredes. fc Klessenl 



20031 see also the reviews by Mac Low & Klessen 2004; 
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Another controversy, related to the control of the 
SFE, refers to the nature of the motions originating the 
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linewidths observed in GMCs and their substructure. 
The latter were initially p roposed to correspond t o 
gravitational contraction by iGoldreich fe Kwanl (|1974l) . 
but this suggestion was qu ickly deemed untenable by 
iZuckerman fc Palmer! (|1974l ). who noted that it would 
imply total Galactic SFRs of the order of the total molec- 
ular gas mass in the Galaxy (~ 10 9 M Q ) divided by the 
typical free-fall time for a GMC (~ 4 Myr), or ~ 250M Q 
yr _1 , an estimate roughly two order s of magnitude larger 
than the observed Galactic SFR. IZuckerman fc Evans 
(|1974[ ) then suggested that the observed linewidths could 
correspond instead to random, small-scal^l turbulent 
motions, a notion that has prevailed until the present. 
However, a number of workers have recently advocated 
a return to the gravitational contraction picture, noting 
that various observational properties of clouds and 
clumps can be well matched by models dominated by 
gravi tational contraction (e.g.. iHartmann fc Burkertl 

21)1)7: iPeretto. rlennebelle fc Andrei 120071: 

Vazquez-Semadeni et al.1 I2009D . Moreover, the no- 



tion of completely random, small-scale turbulent 
motions appears difficult to reconcile with the recent 
realization that the principal component of the velocity 
differences within clouds and clumps at all scales appears 
to be "dipolar", indicative of coh erent motions at the 
scale of the whole cloud or clump (jHever fc Bruntll2"007l ; 
iBrunt et al.ll2009l) . In general, several studies comparing 
simulations and observations have concluded that the 
motions in molecular clouds are consistent with large-, 
rather than smal l -scale driving (lOssenkopf fc Mac Lowl 
l200i IBruntl[200l IPadoan et al.ll2009D . 

If a return to the collapsing scenario is to be 
considered, it is necessar y to somehow avoid the 
IZuckerman fc Palmer! (|1974f ) criticism of it. This is ac- 
tually not so difficult because that criticism neglects 
the internal structure of the GMCs. Recent numeri- 
cal studies of cloud formation by converging streams of 
warm neutral gas in the interstellar medium show that 
the clouds are born turbulent, due to one or more of 
the thermal, thin-shell, and Kelvin-Helmholz instabilities 
dHennebelle fc Peraultlll999HKovama fc Inutsuka 1 1 2000 1 : 
Kovama fc Inutsukal 120021: lAudit fc Hennebelld bpOa 
Heitsch et al.ll2005l : I Vazquez-Semadeni et al J 1 20061) . The 



turbulence is subsonic with respect to the warm gas, 
but supersonic with respect to the cold phase, imply- 
ing large density fluctuations in the latter. In fact, it has 
recently been proposed that molecular clouds may actu- 
ally contain a warmer, atomic s ubstrate in which colder 
mole cular clumps are embedded (jHennebelle fc Inutsukal 
2006). In either case, the molecular cloud contains large, 
nonlinear density enhancements in which the local free- 
fall time is significantly shorter than the cloud's average. 
Thus, once the global collapse begins, the local clumps 
may complete their collapses earlier than the bulk of the 
cloud. They can thus begin forming stars that can begin 
their feedback action on the GMC before it completes 
the bulk collapse. 

In this paper, we present numerical simulations 
aimed at investigating this scenario, in which we use 
the same cloud-formation setup of previous papers 

2 IZuck erman & Evans ( 193) referred to these motions as "lo- 
cal" , and explicitly discarded large-scale coherent motions such as 
gravitational contraction. 



(|Vazquez-Semadeni et~aH 120071: iBaneriee et~aTI 120091) . 
but including a prescription for stellar feedback mimick- 
ing ionization heating from massive stars. With this tool, 
we investigate the effect of the feedback on the global 
SFE of the evolving GMC, as well as the nature of the 
motions in the cloud, in a first effort to shed light on 
these issues. As we shall see, it turns out that the phys- 
ical conditions in the clouds differ significantly from the 
"normal" picture, since accretion of gas from the warm 
diffuse medium is an integral part of the clouds' dynamics 
and evolution, and thus the clouds cannot be considered 
as isolated. 

The plan of the paper is as follows. In ^ we describe 
the numerical code, and the implementation of our star 
formation and stellar feedback prescriptions. In fj3] we 
describe the simulations, and in $4]we describe the results 
concerning the control of the SFE by stellar feedback and 
the nature of the "clouds" themselves. Finally, in <j6]we 
present a summary and some conclusions. 

2. THE NUMERICAL MODEL 
2.1. Heating and cooling 

The numerical simulations used in this work were per- 
formed using the hydrody namics + N-body Adaptive Re- 
finem ent Tree code ART (|Kravtsov et al.Hl997HKravtsovi 
2003). Among the physical processes implemented in 
ART, relevant for our physical problem, are the radia- 
tive heating and cooling of the gas, its conversion into 
stars, ionization-like heating from stellar feedback, and 
self-gravity, from both the stars and the gas. 

We use heating (r) and cooling (A) functions of the 
form 



T = 2.0 x 10~ 26 erg s" 1 



A(T) 



: 10 7 exp 



-1.184 x 10 £ 
T + 1000 



+ 1.4 x 10"VTexp 



-92 



(1) 



(2) 



These functions are fits t o the various heating and coo ling 
processes considered by IKovama fc Inutsuka I (120001) , as 
given by equation (4) of IKovama fe Inutsukal (|2002l ). As 
noted in lVazauez-Semadeni et all ( 2007L h ereafter Paper 
I), eq. (4) in IKovama fc Inutsukal (|2002t ) contains two 
typographical errors. The form used here incorporates 
the necesary corrections, kindly provided by H. Koyama 
(2007, private communication). With these heating and 
cooling laws, the gas is thermally unstable in the density 
range 1 < n < 10 cm~ 3 (cf. Paper I). 

2.2. Star formation and stellar feedback prescriptions 

In ART, star formation is modeled as taking place in 
the coldest and densest regions, defined by T < T$f and 
n > rtgp, where T and n are the local temperature and 
number density of the gas, respectively, and usf and Tsf 
are respectively a density and a temperature threshold. 
We set Tsf = 9000 K, which is easily satisfied by all 
cells with density usf, so in practice our SF condition 
depends on density only. 

A stellar particle of mass is placed in a grid cell 
where these conditions are simultaneously satisfied, and 
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this mass is removed from the gas mass in the cell. There- 
after, the particle is treated as non-collisional, and fol- 
lows N-body dynamics. No other criteria are imposed. 
In each gas cell that satisfies the above criteria a stel- 
lar particle is formed with a mass equal to 50% of the 
gas mass contained in the cell. Since the stellar particles 
are more massive than a single star, each stellar parti- 
cle should be considered as a small cluster, within which 
the individual stellar masses are distributed according to 
some stellar initial mass function (IMF). 

Stellar particles inject thermal energy at a rate E erg 
Myr^ 1 per star with mass greater than 8 M ffl contained 
in the stellar particle. We assume a Mille r fc ScaTol 
(1979) IMF, implying that each stellar particle of mass 
133 Mq produces one 8-Mq star. The energy is de- 
posited in the cell in which the stellar particle is instanta- 
neously located, over a typical OB stellar lifetime, which 
we assume to be 10 Myr. 

It is important to note that, although initially we ex- 
perimented with realistic values of E based on the Ly- 
man continuum fluxes of stars with ma sses between 10 
and 20 M Q (e.g.. lDiaz-Miller et al.lll998l ). we found that, 
because all the energy is deposited in a single cell, and 
the neighboring cells are heated exclusively by conduc- 
tion, rather than by radiative heating, the resulting HII 
regions were not so realistic. Thus, we opted instead for 
taking £ as a free parameter, and adjusting it until we 
obtained realistic HII regions, with temperatures ~ 10 4 
K, diameters of a few parsecs, and expansion velocities 
of a few tens of km s . 

Note also that we resort to the common strategy of 
turning off the cooling in the cell where a stellar par- 
ticle is located, so that the cell can reach realistically 
high temperatures. Otherwise, the cooling can dissipate 
most of the thermal energy deposited in very dense cells. 
In the real ISM this does not occur because the stellar 
heating is applied through ionization, so that the tem- 
perature reached in the star's immediate environment is 
independent of the medium's local density. Instead, in 
the simulations, the cooling depends on the density, and 
the temperature resulting from the balance between the 
stellar heating and the cooling does depend on the den- 
sity. This problem is avoided by turning off the cooling 
in the cell where the stellar particle is located. Note 
that this contradicts claims that the need to turn off the 
cooling ca n be alleviated simply by increasing the resolu- 
tion (e.g.. iCeverino fe~K lvpin 2003). We argue that this 
problem can only be alleviated by performing radiation- 
hydrodynamics simulations. In their absence, we con- 
sider that turning off the cooling is actually a better 
model of the effect of stellar feedback, because it allows 
mimicking the fact that the gas temperature in the HII 
regions is independent of the local density. 

Finally, in order to further constrain the physical con- 
ditions in the HII regions, we also impose a "ceiling" to 
the temperature in the cell containing the stellar particle 
because otherwise, with the cooling off, the temperature 
in the cell might diverge. We set this "ceiling" to 10 6 K. 

Although this procedure is mostly one of trial-and- 
error, we consider it to be the most adequate one for our 
purposes, since it is the HII regions that drive the tur- 
bulent motions in the dense gas in our simulations, and 
so it is them that must have realistic properties, even at 



the expense of a somewhat ad-hoc SF prescription. We 
show a typical HII region in our simulation in Fig. [TJ 

2.3. Refinement 

The numerical box is initially covered by a grid of 128 3 
(zeroth level) cells. The mesh is subsequently refined as 
the matter distribution evolves. The maximum allowed 
refinement level was set to four, so that high-density re- 
gions have an effective resolution of 2048 3 cells, with a 
minimum cell size of 0.125 pc. Cells are refined (halved 
in linear size) if the gas mass within the cell is greater 
than 0.32 M©. That is, the cell is refined by a factor of 2 
when the density increases by a factor of 8, so that, while 
refinement is active, the grid cell size Ax scales with den- 
sity n as Ax oc n -1 / 3 . Once the maximum refinement 
level is reached, no further refinement is performed, and 
the cell's mass can reach much larger values. In par- 
ticular, a stellar particle is formed when a fourth-level 
cell reaches a density hsf = 4 x 10 6 cm~ 3 , or a mass of 
243. 5Mq, again assuming /j, = 1.27 (we use this value be- 
cause we do not follow the actual chemistry, and thus we 
assume the entire box to be filled with atomic hydrogen.) 
Thus, a stellar particle typically has a mass > 122M©. 

Note that, because we use only four levels of re- 
finement, the largest densities arising in the simula- 
tion are by far not sufficiently resolved according to the 
"Jeans criterion" proposed by iTruelove et alj (|1997l ) for 
adaptive- mesh codes. Specifically, at our stellar-particle 
formation threshold density of 4 x 10 6 cm -3 , and assum- 
ing a gas temperature of T ~ 15 K at that density, we 
find that the Jeans length (using the adiabatic sound 
speed) is ~ 0.031 pc, while the minimum cell size, at 
0.125 pc, is roughly 4 times larger. Thus, according to 
those authors, one should expect artificial fragmentation 
to occur in our simulations. However, we do not con- 
sider this to be a problem because we are not concerned 
here with the numbers and masses of the stellar parti- 
cles formed in the simulation, but simply with the total 
amount of mass that goes into stars. 

3. THE SIMULATIONS 

We consider four simulations using the same setup as 
in Paper I, which represents the evolution of a region of 
256 pc per side, initially filled with warm gas at a uniform 
density of no = 1 cm and a temperature Tq = 5000 K, 
implying an adiabatic sound speed c s = 7.4 km s _1 (as- 
suming a mean particle mass fx = 1.27). The whole nu- 
merical box thus contains 5.25 x 10 5 M©. In this medium, 
we set up two streams moving with the same speed 
"inf = 5.9 km s _1 (corresponding to a Mach number of 
0.8 with respect to the unperturbed medium) in opposite 
senses along the rr-direction. The streams have a radius 
of 32 pc and a length of 112 pc each, so that the total 
mass in the two inflows is 2.25 x 10 4 M©. The flows col- 
lide head on at the box's center (see Fig. 1 of Paper I). 
To the inflow velocity field we superpose a field of initial 
low-amplitude turbulent velocity fluctuations, in order to 
trigger the instabilities in the compressed layer. We cre- 
ate this initial velocity fluctuat ion field with a new ver- 
sion o f the spectral code used i n lVazauez-Semadeni et alj 
(1995) and lPassot et alj (|1995f) . modified to run in paral- 
lel in shared-memory architectures. The simulations are 
evolved for about 40 Myr. 
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Fig. 1. — Cross sections of the density (left panel), temperature (middle panel) and y-speed (right panel), shown on the x-z plane, of a 
typical isolated HII region. The scale bar near the top indicates length in parsecs. 



TABLE 1 

Run parameters 



Run 


Wrms 


Feedback 


name 


[km s" 1 ] 




LAFO 


1.7 


off 


LAF1 


1.7 


on 


SAFO 


0.1 


off 


SAF1 


0.1 


on 



The collision nonlinearly triggers a transition to 
the co ld phase, forming a turbulent, cold, dense 
cloud (lAudit fc Hennebelld 120051: IHeitsch et all 120051: 
iHeitsch et alll2006r iVazquez-Semadeni et al.l 120061. con- 
sisting of a complex network of sheets, filaments, and 
clumps of cold gas embedded in a warm diffuse substrate 
(lAudit fc Hennebelldl200aiHennebelle fc Inutsukall2006l: 
iHennebelle fc Auditl I2007TJ . The largest cold structures 
may become gravitationally unstable and begin to col- 
lapse. Eventually, they proceed to forming stars, which 
then heat their environment, forming expanding "HII re- 
gions" that tend to disperse the clumps. 

In the simulations reported here, we vary only two 
parameters: the amplitude of the initial velocity fluc- 
tuations and whether the stellar feedback is on or off. 
We consider a "large-amplitude" (LA) and a "small- 
amplitude" (SA) case for the initial velocity fluctuations, 
for which the three-dimensional velocity dispersions are 
D rms ~ 1.7 km s _1 and u rms ~ 0.1 km s -1 , respectively. 
We thus employ a mnemonic nomenclature for the runs 
using the acronyms LA or SA, followed by F0 or Fl, in- 
dicating that feedback is off or on, respectively. Table [T] 
summarizes the runs considered in the paper. 



4. RESULTS 



4.1. 



Evolution of the simulations 

The simulations performed here behave very similarly 
to previous simulations with similar setups, performed 
with other codes. In particular, our SA runs are very 
similar to run L256Ai>0.1 7 in P aper I and the run pre- 
sented by Bancric e' et al.l ((2009). The main feature of 
these runs is that, because the initial fluctuations are 
very mild, the flow collision creates a large, coherent 
"pancake" of cold, dense gas, which is able to undergo 
gravitational collapse as a whole. This results in the 
formation of a dense, massive, and turbulent region at 
the site where the global collapse finally converges, with 
physical propert ies similar to those of high-m ass star 
forming regions (Vazquez-Sem adeni et al.ll2009f) . How- 
ever, a recent study varying the parameters of the flow 




Fig. 2. — View in projection of the whole numerical box for sim- 
ulation LAF1 at t ~ 31.64 Myr. The box size is 256 pc. The green 
dots indicate stellar particles. The light yellow spots are transient 
dense cores, highlighted by saturation of the color table. In the 
electronic version, this figure shows an animation of the entire evo- 
lution of the simulation up to t = 40 Myr. Records are spaced by 
a time interval At rj 0.14 Myr. 



collision (Ros as- Guevara et al.l I2010h shows that the co- 
herence of the collapse may be lost in the presence of 
stronger initial velocity fluctuations, and the SFE is de- 
creased. In such cases, smaller clouds appeared to be 
less strongly gravitationally bound, with the effect of de- 
creasing the SFE. This feature also happens in our LA 
simulations, in which the cloud formed by the initial flow 
was much more fragmented and scattered over the simu- 
lation volume. As a result, SF also occurs in a much more 
scattered manner, and the SFEs are in general smaller 
in the LA runs than in their SA counterparts. 

However, in general a common pattern is followed by 
all simulations: the transonic converging flows in the dif- 
fuse gas induce a phase transition to the cold phase of 
the atomic gas, which is highly prone to gravitational 
instability. This can be seen as follows. The thermal 
pressure at our initial conditions is 5000 K cm -3 . From 
Fig. 2 of Paper I, it can be seen that the thermal bal- 
ance conditions of the cold medium at that pressure are 
n ~ 130 cm~ 3 , T ~ 40 K. At these values, the Jeans 
length and mass are ~ 7 pc and ~ 64OAf , respectively. 
These sizes and masses are easily achievable by a large 
fraction of the cold gas structures, which can then pro- 
ceed to gravitational collapse and form stars. Moreover, 
the ensemble of these clumps may also be gravitationally 
unstable as a whole, the likelihood of this being larger 
for greater coherence of the large-scale pattern. 

Regions of active star formation form in both sets of 
simulations by the gravitational merging of pre-existing 
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Fig. 3. — Cross-section view through the Central Cloud in run 
SAF1 at t 33 Myr, at which time it has grown to a mass of ~ 
3 X 10 4 Mq . The plane of the image is shown from an inclined line of 
sight for better perspective. The cloud is seen to contain numerous 
HII regions mixed with dense regions. In the electronic version, this 
figure shows an animation of the build-up of this cloud, illustrating 
how it forms by the continued accretion of infalling material from 
the globally collapsing CMC. In the animation, note that many 
of the infalling clumps form stellar particles before reaching the 
center, and are disrupted by the local stellar feedback. However, 
once the Central Cloud is fully assembled, it resists dispersal, and 
forms stars at a high rate. 

smaller-scale clumps, which, altogether, form a larger- 
scale GMC. Figure [2] shows a whole-box image of the 
density field in run LAF1, in projection. In the elec- 
tronic version of the paper, this figure shows an anima- 
tion of this run from t — to t w 40 Myr, illustrating 
the entire evolution of the simulation, from the assembly 
of the cloud, to its advanced star-forming epochs. In the 
animation, subsequent "records" are separated by time 
intervals At 0.14 Myr. 

In the SA runs, the largest star-forming region forms in 
the center of the simulation, due to the coherent collapse 
of the entire sheet-like cl oud formed by the collision. This 
was the region shown in lVazquez-Semadeni et all (2009) 
to exhibit physical conditions typical of actual high-mass 
star-forming regions. We refer to this region as "the Cen- 
tral Cloud" . Figure [3] shows a view of this region at 
t « 33 Myr, a time at which the central cloud has grown 
to a mass of nearly 3 x 10 4 M Q in run SAF1 (cf. Fig. [5]). 
In the case of the LA runs, because star formation occurs 
in a much more scattered fashion, we study two of the 
regions exhibiting the strongest star formation activity, 
neither of which is located at the center of the numerical 
box. These are shown in Fig. |4l and we refer to them as 
Cloud 1 and Cloud 2. 

4.2. Effect of stellar feedback on the SFE and on the 
clouds' evolution 

Our main interest in this contribution is the effect of 
the feedback on the efficiency of the star formation pro- 
cess, and the identification of the mechanism through 
which this effect is accomplished. Figure [5] shows the 
evolution of the dense gas mass and the stellar mass in 
these simulations, with [right panels) and without (left 
panels) feedback. The solid lines refer to the total masses 
in the computational box, while the dotted lines refer to 
the masses in the Central Cloud. Figures [S] and [7] show 
the corresponding plots for Cloud 1 and Cloud 2. Here, 
the solid lines represent the masses for the full simulation 




Fig. 4. — Cross-section view through Clouds 1 (top panel) and 
2 (bottom panel) at t ~ 35 Myr in simulation LAF1. The plane 
of view is located at x = 100 pc for Cloud 1 and at x = 150 
pc for Cloud 2. The dots represent the stellar particles. In the 
electronic edition, these figures show animations of the evolution 
of both clouds from t = 23 to t = 40 Myr. 
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Fig. 5. — Evolution of the dense gas mass and the stellar mass 
for the SA simulations, without (left panels, run SAFO) and with 
(right panels, run SAF1) feedback. The solid lines refer to the total 
masses in the computational box, while the dotted lines refer to the 
masses in a cylinder with length and diameter of 10 parsecs with its 
axis along the ^-direction, and which contains the Central Cloud. 

box, and the dotted, short-dashed, and long-dashed lines 
represent the masses contained in cylinders of length and 
diameter 10, 20, and 30 pc, respectively, enclosing the 
clouds. We do this because the clouds have very compli- 
cated morphologies, with filaments that extend out over 
tens of parsecs and connecting with other clouds (Fig. 
[4]), thus making it virtually impossible to fully enclose 
the "clouds" in any given cylindrical volume. 

It is seen from Figs. [SE that the inclusion of feedback 
(right panels) causes the dense gas mass to be larger and 
the stellar mass to be smaller than in the case without 
feedback in general, even though the total cloud mass 
(dense gas + stars) in the simulations is nearly the same 
in both the cases with and without feedback (Fig. [8]). As 
a consequence, the instantaneous SFE, defined as 



SFE(i) 



M„(t) 



(3) 



where M& c 



M densc (t) + M*(t)' 
2 (t) is the mass of the dense (n > 50 cm~ 3 ) 
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Fig. 6. — Evolution of the dense gas mass and the stellar mass 
for Cloud 1 in the LA simulations, without (left panels) and with 
(right panels) feedback. The solid lines refer to the total masses 
in the computational box. The other lines refer to the masses in 
cylinders of length and diameter 10 pc (dotted lines), 20 pc (short- 
dashed lines), and 30 pc (long- dashed lines). 

gas in the simulation, and M*(i) is the stellar mass, nat- 
urally decreases upon the inclusion of feedback in both 
sets of simulations (Fig. |9j). Note that in eq. J3]) we have 
explicitly written out the time dependence of the quan- 
tities involved. 

The SFE is seen to be reduced by a larger factor (~ 
10 x) in the case of the LA runs, in which the collapse is 
less focused and less massive, than in the case of the SA 
ones (~ 3x), in which the opposite is true. In addition, 
in Fig. [10] we show the evolution of the SFE at the level of 
the clouds. The left panel shows the evolution of the SFE 
for the Central Cloud in the SA runs. The middle panels 
show the corresponding plots for Cloud 1 and Cloud 2 in 
the LAFO run (without feedback), and the right panels 
show the SFEs in the LAF1 run (with feedback). Again 
we see a trend for the less massive cloud (Cloud 1) to 
suffer a greater reduction of its SFE (by a factor of ~ 20, 
from ~ 60-90% to ~ 3-4%) than the more massive one 
(Cloud 2, by a factor of - 10, from ~ 70-80% to ~ 7-8%). 
We discuss the possible cau ses for this mass-dependent 
effect of the feedback in Sec. 14.3.21 

The factor by which the SFE is reduced upon the inclu- 
sion of the feedback in the simulations is plotted versus 
the system's mass in Fig. [TTJ both for the full amount of 
dense gas in the simulations, and for each of the clouds we 
have been considering. We see that two sets of points are 
clearly defined in this plot, one for the clouds, and one 
for the simulations. In both cases, however, the trend 
of larger reduction factor at smaller dense gas mass is 
clearly observed, although, at the level of simulations, 
we see that their masses are not very different. Thus, 
in this case the different reduction factors must include 
a contribution from the larger degree of fragmentation 
occurring in the LA simulations due to the larger am- 



FlG. 7. — Evolution of the dense gas mass and the stellar mass 
for Cloud 2 in the LA simulations, without (left panels) and with 
(right panels) feedback. The solid lines refer to the total masses 
in the computational box. The other lines refer to the masses in 
cylinders of length and diameter 10 pc (dotted lines), 20 pc (short- 
dashed lines), and 30 pc (long- dashed lines). 

plitude of the initial turbulent fluctuations, beginfigure 

In order to compare the SFEs in our simulated clouds 
with those of real molecular clouds, it is convenient to 
note tha t regions of low-mass star formation, such as 
Taurus ( Goldsmit h et al.] |2008) or the Chamaeleon II 
dark cloud (jSpezzi et al.l 12008! ) generally have low SFEs 
(~ 1- 5%), while cluster forming cores have SFE — 30- 
50% (ILada fc Ladal[200l . Thus, we can check whether 
the SFEs and SFRs of our three clouds follow the same 
trend. Figure [T2l shows the evolution of the SFRs for our 
three clouds. From these plots, we find that the average 
SFR of the Central Cloud starting from t = 32 Myr (the 
time at which a large, roughly stationary SFR sets in) is 
(SFR) - 1450 M Q Myr" 1 , while for Clouds 1 and 2, we 
find (SFR) — 30 and 60 M Myr -1 during their entire 
star-forming stages, respectively. From here, and using 
the instantaneous SFE at t = 40 Myr from Fig. QJJ1 we 
can then plot the SFE versus the SFR. This is shown 
in Fig. [U in which both the SFR and the SFE have 
been multiplied by an extra factor of 0.5, to respresent 
the fact that the efficiency within the stellar particles, 
which are created at a density of n = 4 x 10 6 cm~ 3 in 
our simulations, is still smaller than unity. We take 0.5 
as a representative value. 

We see that the Central Cloud has values of the 
SFR and the S FE comparab l e to those of cluster- 
forming cores (ILada fc Ladal [2003). Specifically, 
iVazquez-Semadeni et al.l (|2009f) estimated an SFR > 
25OM Myr" 1 for the Orion Nebula Cluster (ONC). 
This calculation u sed an estimated age s pread of the 
ONC of < 2 Myr ( Hil lenbrand! [I997h . and iTobin et~aTl 
(2009) 's result of there being 1613 stars in the ONC. 
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Fig. 8. — Evolution of the total (dense gas + stars) mass for the 
four simulations. The SA simulations are shown in the top panel, 
while the LA runs are shown in the bottom panel. The black, solid 
lines refer to simulations with feedback and the red, dotted lines 
represent runs without it. The colors are shown in the electronic 
version only. 
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Fig. 9. — Evolution of the instantaneous star formation efficiency 
(SFE), as defined in eq. J3J, in the full simulation box in the four 
runs. 



Taking this number as a proxy for the total stellar 
producti on of this region, and a mean stellar mass of 
0.3 M Q ([Hillenbrand fc Carpenter! I2000D . this implied a 
total stellar mass of ~ 500Af Q . On the other hand, 
iKrumholz fc Tanl (l2007fl quote a total stellar mass i n the 
ONC of ~ 4600M© (Hillenb rand fc Hartmanr] ll99"l and 
an age spread of ~ 3 Mvr (|Tan et al.ll2006l ). implying an 
SFR - 153OAf Myr- 1 . These values bracket the SFR 
we measured for our Central Cloud. 

Concerning the SFE, Fig. [T3] shows that the Central 
Cloud has SFE ~ 10%, which is c omparable to that of 
the Orion A cloud (|Carpenterll2000D . in which the OMC- 
1 clump is contained. Thus, our Central Cloud may be 
compare d to the Orion A cloud, and its dense core, dis- 
cussed in lVazquez-Semadeni et all (|2009). is comparable 
to the OMC-1 clump. 

On the other hand, Clouds 1 and 2 are seen in Fig.[T3lto 
have SFRs comparable to those of low-mass star-forming 
clouds, ~ 10-100 M Q Myr" 1 , such as Taurus, Cha II or 
Lupus, in which th e SFE is known to be lower, ~ 1-4% 
(jSpezzi et al.l 12008). This compares very well with the 
SFEs we report for our Clouds 1 and 2 in the same figure. 
We conclude that Clouds 1 and 2 are good models of 
low-mass star-forming regions, while the Central Cloud 
is a good model of a ma s sive S F region, as discussed in 
iVazquez-Semadeni et al.1 (|2009l ). 

4.3. The physical processes acting on the clouds 
4.3.1. Cloud "destruction" 

Up to the 40-Myr time to which we have evolved our 
simulations, the three large clouds (the Central Cloud in 
run SAF1 and Clouds 1 and 2 in run LAF1) do not show 
any instances of the dense gas mass reversing its increas- 
ing trend and beginning to decrease due to the feedback, 
as can be seen in Figs. [5] through [71 Apparently, the 
HII region-like feedback we use is unable to overwhelm 
the large gravitational potential well of these clouds and 
their enveloping "atomic" gas reservoirs. 

Instead, complete destruction seems to be able to occur 
in small clumps. This can be seen in the animation cor- 
responding to Fig. [3j in which various small clumps are 
seen to be destroyed by their stellar products. Three par- 
ticularly conspicuous ones are, first, the one that forms 
a stellar particle at the very starting frame (record 179) 
of the animation, slightly above and to the right of the 
screen's center. Next, another stellar particle appears in 
the clump almost at the left border of the frame, about 
2/3 of the way from bottom to top, at record 187. Fi- 
nally, a third particle appears at record 189, slightly be- 
low and to the right of the screen's center, as the result of 
the collision of two clumps. In all of these cases, a single 
stellar particle is formed (~ 120M©), and the clump is 
destroyed. It is worth noting that actually, the expansion 
of the HII regions formed produces new clumps from the 
material collected around it, but these new clumps cither 
disperse, or simply do not form new stars. 

Thus, we conclud e, similarly to 

IKrumholz. Matzner. fc McKed (|2006f ). that small 
clouds ("clumps") are rapidly destroyed, while large 
clouds may survive for longer times. However, our 
clouds exhibit a fundamental difference with respect to 
the model considered by those authors, namely that the 
clouds are accreting in general. In the next section we 
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Fig. 10. — Evolution of the instantaneous star formation efficiency (SFE) in the dense clouds in the simulations. The left panel shows 
the SFE for the Central Cloud of the SA runs, with and without feedback. The measurements refer to a cylindrical box with a diameter 
and a length of 10 pc. The middle panels show the SFE for Clouds 1 and 2 in the LAFO simulation (without feedback), for three different 
cylindrical boxes, of length and diameter indicated in the labels. The gaps in the curves for the smaller cylindrical boxes correspond to 
times when the stellar particles migrate out of them, and no new particles have been formed. The right panels show the SFE for Clouds 1 
and 2 in the LAF1 simulation (with feedback). 
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Fig. 11. — Reduction factor of the SFE at 40 Myr upon inclu- 
sion of the feedback in the simulations, both for the three clouds 
as well as for the entire mass of dense gas in the simulations. The 
error bars in the masses indicate the range of values taken b y th e 
systems once their initial assembly has finished (see Figs. I5l71l . 
The error bars in the reduction factor for Clouds 1 and 2 denote 
the maximum and minimum values of the ratio SFE(w/o feed- 
back)/SFE(w/feedback), using the SFE data for the three cylinder 
sizes of 10, 20 and 30 pc. 

now discuss this feature. 

4.3.2. Accretion vs. feedback 

One crucial feature in all our simulations is that the 
clouds are accreting material from the surrounding dif- 
fuse medium. This is fundamentally different from mod- 
els in which the clouds are isolated entities, in rough bal- 
ance between their self-gravity and the turbulent pres- 
sure, possibly driven by the stellar feedback. The accre- 
tion competes with SF and stellar feedback in regulating 
the cloud's mass and coherence, with important conse- 
quences. First of all, this implies that simple observa- 
tional estimates of the SFE in GMCs based on measur- 
ing the stellar mass and dividing it by the cloud's mass 
may be failing to take into account the additional "raw 
material" for SF contained in the part (or the whole) of 
the atomic envelope of the clouds that will eventually be 
incorporated into the GMC. 

Second, the competition between feedback and accre- 
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Fig. 12. — Evolution of the stellar particle formation rate, aver- 
aged over 2-Myr intervals, of the Central Cloud (top panel), Cloud 
1 (middle panel) and Cloud 2 (bottom panel). For this plot we 
only use the 10-pc cylindrical volumes for Clouds 1 and 2. Because 
our stellar particles form at a density of n = 4 X 10 6 cm -3 , actual 
star formation rates should be a factor of 2-3 lower. Gaps indicate 
periods over which no new stars are formed. 
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Fig. 13. — SFE versus stellar-particle formation r ate for the three 
clouds in the simulations, using the data from Figs.lTolandlT2l The 
values of both the SFR and the SFE have been multiplied by a 
factor of 0.5, representative of the still-lower-than-unity efficiency 
within our stellar particles. 

tion may explain our observation from Sec. 14.21 that 
cases with feedback are characterized by larger dense 
gas masses and smaller stellar masses than their coun- 
terparts without feedback. The smaller stellar mass is 
not surprising, as the obvious effect of stellar feedback 
is to reheat the cold, collapsing star-forming gas, thus 
reducing the SFR. However, the larger cold gas mass in 
the presence of feedback is indeed surprising, since both 
gas consumption by SF and the "ionization" by stellar 
feedback act to reduce the dense gas mass. Our result 
implies that the rate of dense gas consumption by SF 
far outweighs its rate of destruction by stellar feedback, 
so that the net effect of reducing the SFR is to allow a 
larger amount of dense gas to be collected by the accre- 
tion. This scenario is supported by Fig. [5J which shows 
the total (dense gas + stars) in the clouds in the two 
sets of simulations. It is seen that the total cloud mass is 
nearly the same with and without feedback, suggesting 
that the total cloud mass is mainly determined by the ac- 
cretion, while the ratio of dense gas to stellar mass seems 
to be mainly determined by the feedback. 

These results moreo ver support the suggesti o ns by 
lElmegreenl (|2007l ) and IVazquez-Semadeni et ail (|2009f ) 
that the majority of the gas in a GMC is not partici- 
pating of the SF process at any given time. The latter 
authors suggested that this is how local regions of SF 
may have a very high specific SFR (~ (10 Myr)" 1 ), while 
that of their whol e parent GMCs may be m uch smaller 
(~ (300 Myr)" 1 ; IMcKee fe Williamsl[l99l . Because it 
is injected by the newly formed stars, the stellar feedback 
acts preferentially on gas that is about to form stars next. 
This allows an efficient suppression of SF, through only 
a modest fraction of the total available dense gas being 
destroyed. That is, if SF is a highly localized process, 
then it is possible to achieve significant reductions of the 
SFR with only a modest reduction of the total amount of 
dense gas, by targeting the destruction precisely at the 
star-forming gas. 

This mechanism may also explain the trend observed 
in Sec. 14.21 that the SFE is more strongly reduced by 
the feedback in cases where the collapsing gas mass is 
smaller. This may be understood as a consequence of the 



fact that stellar feedback is localized, while the accretion 
is extended, and more so for greater mass of the globally 
gravitationally unstable region that will form the cloud. 
Thus, as we observe in the animation corresponding to 
Fig. |3l the low-mass fragments that are undergoing SF 
while on their way to the central site of the global collapse 
in run SAF1 are easily destroyed by their stellar activity. 
Instead, the massive Central Cloud is not destroyed, as it 
continues to accrete mass from large distances at a pace 
that outweighs the local destruction by stellar feedback. 
Clouds 1 and 2 in run LAF1, which do not involve such 
an extended collapse, are intermediate cases, in which 
the clouds are not destroyed by the feedback, but the 
latter is more efficient in reducing the SFE. 

Thus, we can conclude that the efficacy of the feed- 
back in destroying the cloud is maximal when its region 
of influence is comparable to the spatial extent of the 
infall, and is progressively reduced as the latter involves 
progressively larger coherence lengths. 

4.4. Physical conditions in the clouds 

Here we consider the global physical properties of the 
clouds. We postpone a discussion of the properties of 
individual clumps within the clouds for a future study, 
to be performed at higher resolution. One such study, 
including magnetic fields, although not i ncluding stellar 
feedba ck, has recently been presented by iBaneriee et al.l 

4.4.1. Density PDFs 

It is important to determine the physical conditions 
in our clouds in order to assess their degree of realism. 
One basic diagnostic is the probability density function 
(PDF) of the density field. Although it is well established 
that in iso thermal flows the density PDF takes a lognor- 
mal form (lVazquez-Semadeni|[r9 94; Pad oan et al.l 119971; 
i Passot fe Vazquez-Semadenilll998t iNordlund fe Padoanl 
Il999t lOstriker et al.lll999L 120011: iFederrath et al.ll2008f) . 
for non-isothermal flows the expectation is in gen- 
eral different, with a near-power-law tail develop- 
ing a t high densities for flows s o fter th an isother- 
mal (IPassot fc Vazque z-Scmadcni 1998; Scal o et all 
1998t INordlund fc Padoanl Il999t IGazol et all 120051: 
de Av illez fc Brcits chwerdtl 12004 although see, e.g., 
Wada & Norman 2001, 2007), and a bimodal form aris- 
ing for thermally b i stable flows (IVazquez-Semadeni et al.l 
120001 IGazol et al.l 120051 lAudit fc Hennebelld l2005T~ 5ee 
also the review by Vazquez-Semadeni 2009). 

Figure Q3] shows the density PDFs for the whole sim- 
ulation box of the SA runs top panels and of the LA 
runs bottom panels. The right panels show the entire 
density range, while the left panels show the PDF for 
the dense gas (n > 100 cm~ 3 ) only. The whole-range 
PDFs exhibit the bimodality typical of thermally bistable 
flows, although the high-density tail is seen to exhibit 
an excess over the power-law in both the cases with 
and without feedba ck. This is probably due to the ac- 
tion of self-gravity (IKlessenl [20001 IDib fc Burkertl 12001 
IVazquez-Semadeni et al.l 120081 ). In addition, the cases 
with feedback show a slight excess over the cases with- 
out it at densities 10 3 < n < 10 5 cm -3 , probably due 
to the formation of compressed regions by the expand- 
ing "HII regions". The high-density PDFs, being just 
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Fig. 14. — Density PDFs for the entire simulation boxes of the 
SA runs (top panels) and the LA runs (bottom panels). The right 
panels show the PDFs for the entire range of densities, while the left 
panels show the PDFs of the dense gas only (n > 100 cm -3 ). The 
red, dotted lines refer to the simulations without feedback, while 
the black, solid lines show represent the simulations with feedback. 

the n > 100 cm -3 tail of the whole-range PDFs, only- 
show more detail of the very-high density gas, up to the 
star- forming density of 4 x 10 6 cm~ 3 , but with the same 
shape of the curves as in the whole-range PDFs. 

In order to see the PDFs at the locations of the 
three main clouds we have studied (the Central Cloud 
in the SA runs and Clouds 1 and 2 in the LA 
runs), we show in Fig. [15] their corresponding density 
PDFs. It is noteworthy that the PDFs again show a 
roughly power-law shape at high densities over three 
to four orders of magnitude in density, in agreement 
with the expectation for softer-than-isothermal flows 
dPassot fe Vazquez-Semadenilll998l : lNordlund fc Padoanl 
1999). This is at odds with results from numerical sim- 
ulations of turbulent isothermal gas in closed boxes, but 
then again our clouds are not isothermal. Instead, they 
are characterized by an effective polytropic equation of 
state P cx pj s , with 7 c ff varying from 0.8 to nearly unity 
for n > 100 cm" 3 , as can be seen in Fig. [16] Real in- 
terstellar molecul ar gas is not exp ected to be e xactly 
isothermal, either dScalo et alJl99ctlSpaans fc Sflkl l2000l: 
Uappsen et al.ll2005l) . Moreover, our clouds are immersed 
in a diffuse, warme r medi um w hich, as proposed, e.g., 
by iLi fe Goldsmith] (|2003| ) and iHennebelle fc Inutsukal 
(2006), infiltrates the clouds to some degree. Thus, the 
density PDFs for the volumes containing our clouds also 
exhibit the bimodal shape characteristic of thermally 
bistable flows, and have low-density tails that extend 
down to the WNM regime, contrary to what happens 
in isothermal simulations. 

4.4.2. Velocity dispersions and virial masses 

Finally, we investigate the global velocity dispersion 
(A v ) in the Central Cloud and Clouds 1 and 2. This 



Fig. 15. — Density PDFs for the three main clouds in the simu- 
lations: the Central Cloud in the SA runs (left panel), and Clouds 
1 (middle panel) and 2 (right panel) in the LA runs. The PDFs 
for the Central Cloud are computed in a cylinder with length and 
diameter equal to 10 pc, while the PDFs for Clouds 1 and 2 are 
computed in cylinders of length and diameter equal to 10, 20 and 
30 pc. Red lines indicate cases without feedback, and are dis- 
placed downwards by a factor of 10 for better viewing in the cases 
of Clouds 1 and 2. Black lines indicate cases with feedback. 

is shown in Fig. [T7J The solid lines show the density- 
weighted value, while the dotted lines show the volume 
weighted value. Black lines denote cases with feedback, 
and red lines correspond to cases without it. 

For the Central Cloud and Cloud 2, which are the two 
most massive ones, the density-weighted velocity disper- 
sion, which highlights the dense gas, decreases upon the 
inclusion of feedback. Without feedback, in the Cen- 
tral Cloud, it reaches very high values at the end of 
the run (~ 15 km s _1 , corresponding to Mach numbers 
M s ~ 50-75), which in fact are significantly larger than 
typical values for cl oud complexes of comparable mass 
(M_> 10 4 M Q ; e.g., iDame et all 119861 : iRathborne et al.l 
2009). Instead, in the run including feedback, A v reaches 

values ~ 5-6 km s _1 (rms Mach number, M s ~ 20), in 
much better agreement with typically observed values. 
On the other hand, the volume-weighted velocity disper- 
sion, which highlights the less dense gas, is seen in Fig. 
1171 to increase upon the inclusion of feedback. In Clouds 
1 and 2, which are less massive and more scattered than 
the Central Cloud, A v does not achieve exceedingly large 
values in any case. 

These results can be understood as a consequence of 
the fact that the dense gas acquires its largest velocities 
in the case of free-fall collapse. However, the collapse 
flow is dismantled in its final (fastest) stages by the stel- 
lar feedback, so that it is the high-velocity dense gas that 
is preferentially destroyed by the feedback. On the other 
hand, this gas becomes warm, diffuse, high-velocity ex- 
panding gas, which is the one highlighted by the volume 
weighting. These results again reinforce the notion that 
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Fig. 16. — Phase-space diagrams of temperature vs. density for 
the Central Cloud (left panels), Cloud 1 (middle panels) and Cloud 
2 (right panels). The top panels show the cases without feedback, 
and the bottom panels show the cases with feedback. It can be seen 
that, for the majority of the points, the temperature only varies 
between ~ 40 and 10 K for 100 < n < 10 6 cm -3 . The exceptions 
are locations recently heated by new stars. 

the feedback is mostly applied on the gas that is closest 
(but not quite there yet) to forming stars. 

Once we have determinations of the velocity dispersion 
in the cylinders containing the clouds, it is natural to 
measure the virial mass (Af v j r ) of the clouds and compare 
it with their real mass M, in order to check whether they 
look "virialized" . We compute the virial mass through 
the standard formula 
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R\ ( Av 



pcy \ km s 1 



Mr. 



(4) 



(see, e.g., iCaselli et all 1200$ iTachihara etH] l2002t 
Kle ssen et all 120051) . For the real cloud mass we con- 
sider the sum of the gas (dense + diffuse) and stellar 
masses in the volume being considered. It is impor- 
tant to note that, with this prescription, variations in 
the value of the ratio come exclusively from the estimate 
of the virial mass, because it does depend on the weight- 
ing used, while the actual cloud mass is independent of 
the method of weighting. 

Figure IT51 shows the evolution of the ratio M/M v - a for 
the three clouds, showing the cases with and without 
feedback, and using volume- and density-weighting for 
the calculation of A v . In all cases we observe that, at 
advanced stages of the evolution, when the clouds have 
been fully assembled (t > 25 Myr), this ratio is closest to 
unity when feedback is included and density- weighting is 
used in the calculation of A v . The cases with no feedback 
and density weighting are systematically lower than unity 
by factors ~ 2-10. The cases with volume weighting are 
more strongly fluctuating, but it is noteworthy that the 
ratio often takes values close to unity anyway. 
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Fig. 17. — Evolution of the velocity dispersion <r in the three 
clouds. The solid lines show the density-weighted value, while the 
dotted lines show the volume weighted value. Red lines correspond 
to cases without feedback, and black lines correspond to cases with 
feedback. 



From these results, we conclude that, in the case with 
feedback, the clouds appear to be in a pseudo- virialized 
state with respect to the velocity dispersion of the dense 
gas (highlighted by the density weighting). In this 
state, there is an approximate force balance between self- 
gravity and feedback driving. However, it differs from 
true virialization because of the presence of mass sources 
(the accretion from the diffuse environment) and sinks 
(the consumption by SF and the destruction by stellar 
feedback) in the system. These seem to self-regulate, so 
as to be capable of maintaining an approximately con- 
stant cloud mass while accretion persists. On the other 
hand, when the velocity dispersion of the warm gas is 
taken into account (still in the case with feedback), the 
clouds appear to be slightly sub-virial, suggesting that 
the newly formed warm gas in the HII regions is capable 
of escaping the cloud. 

In the cases without feedback, the clouds also appear 
to be sub-virial even with respect to the density- weighted 
virial mass estimate. In this case this appears to be a con- 
sequence of the dense gas being in free-fall collapse, and 
into a potential well produced not only by its own self- 
gravity, but also by t he previously-formed stars that have 
fallen there too (see iVazquez-Semadeni et al.ll2009f ). 
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Fig. 18. — Evolution of the ratio of actual to virial mass for the 
three clouds. Red lines correspond to the runs without feedback. 
Black lines correspond to the runs with feedback. Solid lines cor- 
respond to the density-weighted velocity dispersion, while dotted 
lines correspond to the volume-weighted one. 



5. CAVEATS AND LIMITATIONS 

Of course, our simulations are not free of caveats and 
limitations. Two that stand out are the resolution and 
the nature of the feedback that we have included. Con- 
cerning the resolution, in this paper we have used a rela- 
tively limited one, in order to speed up the simulations, 
at the ex pense of not fu l filling the Jeans criterion pro- 
posed by iTruelove et al.l (|1997t ). However, as explained 
in Sec. 12.31 we do not consider this a problem for our 
study, since we have restricted it to the global properties 
of the clouds, rather than their detailed structure, and we 
have avoided addressing issues related with the fragmen- 
tation of the final stages of collapse in the clouds, such 
as the mass distribution of the cores and stellar prod- 
ucts. We plan to perform a higher-resolution study in a 
forthcoming paper, in which we can address these issues. 

The second limitation is the nature of the stel- 
lar feedback that we have considered, which in this 
study has been restricted to local heating represent- 
ing the ionizing radiation from massive stars , sim i- 
larly to the approach used by, e.g.. lYorke et all (119821) : 
iVazquez-Semadeni et al.l (|1995h : lPassot et aiT i 19951 ). ne- 
glecting oth er sources such as outflows from stars of all 
masses (e.g.. lNorman fe Silk 1980l: lLi fe N akamura1 l200d ; 



Wang et al.ll2009ir We do this in part for the numerical 
simplicity of the approach, and in part because here we 
have been mainly interested in the SFE and and cloud 
evolution at the scale of GMCs, for whic h the expansio n 
of HII regions is likely the main driver ()Matznerl 12002). 
However, the neglect of bipolar outflows may introduce 
a non-negligible bias in our finding that smaller-scale 
clumps are destroyed by the feedback, and this should 
be confirmed by future simulations that can better re- 
solve these objects and include bipolar outflows. 

6. SUMMARY AND CONCLUSIONS 

In this paper we have presented a numerical investiga- 
tion of the evolution of dense ("molecular") clouds, start- 
ing from their formation by transonic compressions in the 
warm neutral medium, followed by a phase transition to 
the cold neutral medium, the onset of gravitational col- 
lapse, star formation, and, finally, energy feedback from 
massive-star ionizing radiation and the formation of ex- 
panding HII regions. 

A crucial difference between the results from 
our simulations and other models and theories of 
the self-regulat ion of sta r form a tion and it s ef- 
ficiency (e.g.. iWhitworthl 119791: lElmegreenl 119831: 



Franco. Shore fe Tenorio-Tagld 1 1991 iMcKee fc Tanl 
2003t iKrumholz fc McKed 120051 ) is that the clouds in 



our simulations are in general accreting material at high 
rates from the surrounding diffuse medium, rather than 
having a fixed mass. This implies that the material 
making up a cloud is constantly being renewed over time 
because, on the one hand it is accreting fresh gas, and on 
the other it is losing mass to SF and stellar ionization. 
The star-forming regions can be considered to be not 
objects, but rather the loci where the gas is just "passing 
through" , from the diffuse-gas state to the "star" state, 
similarly to the nature of a flame, which is the locus of 
the gas undergoing combustion in a candle, with fresh 
air entering at its base, burning while it transits through 
the flame, and the exhaust gases leaving at the top. 

Witin this framework, we have found that the SFE 
in the clouds is readily decreased by feedback to levels 
consistent with observational determinations at all times 
during the clouds ' evolution up to the maximum integra- 
tion time of 40 Myr that we have considered. This is a 
significant improvement over our previous non-magnetic 
studies of the SFE in the context of cloud evolution 
without feedback, in which th e SFE at late times is of- 
ten found to be excessive (e . g.. IVazauez-Semadeni et al.l 
l2007t iRosas-Guevara et al.l |2010TT However, we have 
found that the reduction factor upon the inclusion of 
feedback is an inverse function of the dense gas mass in 
the system, and of the degree of coherence of the global 
collapse, as illustrated in Fig. [TTJ This result is com- 
plemented by the additional observation that low-mass 
clumps are fully destroyed by the feedback. Together 
they imply that at some point, below a certain critical 
collapsing mass, the appearance of a stationary state may 
not occur any more. This is in qualitative agre e ment with 
the results of IKrumholz. Matzner. fc McKed <|2006f l. al- 
though their model does not include mass accretion onto 
the cloud. We plan to present one such model in a future 
paper. 

In general, our results may provide an explanation 
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for the observational fact that regions of massive-star 
formation, which are themselves more massive than re- 
gions of low mass-star f ormation, appear to have higher 
SFEs (La da fc Lad a 2003) t han regions for ming low- and 
intermediate-mass stars (jEvans et al.ll2009I L Essentially, 
our results imply that the more massive the region, the 
less effective the feedback is in reducing the SFE. Specif- 
ically, these results indicate that the role of the feedback 
is not the same in all clouds, but rather depends on the 
initial conditions of the large-scale collapse that produces 
each individual region. Regions in which the amount of 
mass involved in the collapse overwhelms the destruc- 
tive action of the feedback may reach a stationary state 
that appears virialized, but in which there is actually 
a countinuous processing of material, since on the one 
hand the cloud is accreting mass from the outer infalling 
material, and on the other hand it is losing mass through 
consumption by SF and evaporation by the feedback. In 
this case, termination of the SF episode probably re- 
quires the termination of the gas reservoir involved in 
the global collapse, and is independent of the feedback. 
However, this need not be in contradiction with the low 
observed global efficiency of SF, if only a small fraction of 
a CMC's mass is involved in the collapse at any moment 
in time, with the rest being either supported ( by, e.g., the 
magn etic field), or in the process of dispersal (Elmcgreen 
2007). Alternatively, the SF episode may be terminated 
upon the initiation of supernova events, which we have 
not included in the present study. 

An additional consequence of the mechanisms de- 
scribed in this paper is that the most massive clouds 
appear to actually contain more mass when feedback is 
included than when it is not. This suggests that, first, 
unimpeded SF is more efficient at removing mass from 
the dense gas phase to deposit into stars, than the evapo- 



ration of dense gas by the stellar feedback, so that, when 
the latter inhibits SF, the accretion onto the cloud ac- 
cumulates larger amounts of dense gas than otherwise. 
Second, this indicates that the deposition of the feedback 
energy into the clouds is accurately targeted to gas that 
is already on the verge of forming stars. If this gas is a 
minor fraction of the total dense gas in a cloud (i.e, SF is 
a highl y spatially intermittent phenom enon in molecular 
clouds nVazquez-Semadeni et al.M2009h . then its destruc- 
tion is a highly efficient way to reduce the SFE without 
destroying a large fraction of the cloud. 

Finally, we have also investigated the density PDF in 
the volumes containing the clouds in our simulations, 
finding that in general they retain the bimodal shape 
characteristic of thermally bistable flows. The PDF 
of the dense gas exclusively shows no turnover at low 
densities, indicating that low densities are most abun- 
dant. The high-density tails of the PDFs have power-law 
shapes, as expected for softer-than-isothermal flows, al- 
though this last result may be biased by our usage of a 
single cooling law appropriate for atomic gas, and extrap- 
olating it to molecular-gas densities, rather than using a 
specific molecular cooling law. We expect to address this 
shortcoming in future works. 
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